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ABSTRACT 

We present maps of the cosmic large-scale structure around the twelve most distant galaxy 
clusters from the Massive Cluster Survey (MACS) as traced by the projected surface density 
of galaxies on the cluster red sequence. Taken with the Suprime-Cam wide-field camera on 
the Subaru telescope, the images used in this study cover a 27 x 27 arcmin 2 area around each 
cluster, corresponding to 10 x 10 Mpc 2 at the median redshift of z = 0.55 of our sample. We 
directly detect satellite clusters and filaments extending over the full size of our imaging data 
in the majority of the clusters studied, supporting the picture of mass accretion via infall along 
filaments suggested by numerical simulations of the growth of clusters and the evolution of 
large-scale structure. A comparison of the galaxy distribution near the cluster cores with the 
X-ray surface brightness as observed with Chandra reveals, in several cases, significant offsets 
between the gas and galaxy distribution, indicative of ongoing merger events. The respective 
systems are ideally suited for studies of the dynamical properties of gas, galaxies, and dark 
matter. In addition, the large-scale filaments viewed at high contrast in these MACS clusters 
are prime targets for the direct detection and study of the warm-hot intergalactic medium 
(WHIM). 

Key words: galaxies: clusters: general - large-scale structure of Universe - X-rays: galaxies: 
clusters 



1 INTRODUCTION 

As the largest gravitationally bound constituents of matter, galaxy 
clusters provide important clues about the distribution and funda- 
mental properties of galaxies, gas, and dark matter. Studying their 
characteristics over a range of redshifts can yield crucial insights 
into the overall structure and evolution of the cosmos. Theoret- 
ical considerations and the results of numerical simulations (e.g. 
Colberg et al. 2000; Evrard et al. 2002) predict that clusters grow 
through accretion and merger events and reside within filamentary 
structures extending up to 10/ h Mpc in size. While growth through 
cluster mergers is a relatively well studied phenomenon, the exis- 
tence and relevance of filaments for cluster and galaxy evolution 
remains much less well understood, primarily because of observa- 
tional challenges created by the very large scales and low densities 
involved. 

Numerical simulations (e.g. Cen & Ostriker 1999, Dave et 
al. 1999) suggest that a large fraction of the baryons in the local 
universe takes the form of highly diffuse gas at a temperature of 
10 5 — 10 7 i^ - the so-called "Warm-Hot Intergalactic Medium" 
(WHIM). The temperature and diffuse nature of the WHIM makes 
it very difficult to detect directly, the currently most promising tech- 
niques being very deep observations at extreme ultraviolet (EUV) 
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and soft X-ray wavelengths, to detect high-ionization metal absorp- 
tion lines and diffuse bremsstrahlung emission, respectively (e.g. 
Soltan, Freyberg & Hasinger 2005, Nicastro et al. 2005). Since sim- 
ulations also predict that the WHIM is distributed along large-scale 
filaments around galaxy clusters (Dave et al. 2001) and may be de- 
tectable in observations of the Sunyaev-Zel'dovich Effect (Hallman 
et al. 2007), promising targets for the direct detection and subse- 
quent study of the WHIM could be unveiled by a systematic survey 
of the large-scale structure surrounding massive galaxy clusters. 

Large optical surveys of the nearby universe such as 2dF and 
SDSS have allowed detailed investigations of the galaxies within 
the vast cosmic web consisting of clusters, superclusters, voids, 
walls, and filamentary structures (e.g. Peacock et al. 2001; Zan- 
divarez, Merchan & Padilla 2003; Erdogdu et al. 2004; Einasto et 
al. 2007). Because of the limited depth of the underlying imaging 
and spectroscopic data obtained by these surveys, studies of low- 
density environments have focused on structures at relatively low 
redshift (z ~ 0.1 — 0.2). A search for filamentary structure around 
the vertices of the "cosmic web" would, however, be particularly 
rewarding at higher redshifts where the contrast is greatly enhanced 
by the smaller angular size of such features. In addition, observa- 
tions at higher redshift would establish an important baseline for 
evolutionary studies of this environment. Galaxy clusters at higher 
redshifts are expected to be less evolved and reside within the dens- 
est regions of the precursors of the structure seen locally. A large 
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fraction of the matter, dark and luminous, in these structures will 
eventually become part of virialized cluster cores through merging 
and infall along the filaments. Hence a multitude of astrophysical 
topics can be addressed by studying the large-scale surroundings 
of distant clusters, including the impact of environment on star for- 
mation in galaxies, the relative importance of gradual infall along 
filaments vs. merging for both galaxy and cluster evolution, and the 
existence and properties of the WHIM. 

The study described in the following looks for evidence of 
filamentary structure in a complete sample of twelve moderately 
distant, very massive galaxy clusters in order to establish obser- 
vationally the frequency and prevalance of such structure around 
the densest nodes of the cosmic web. By targeting galaxy clusters 
at intermediate redshift (z ~ 0.5), we observe filaments on angu- 
lar scales of typically 10 arcmin, which makes them accessible to 
efficient observational sampling in wavebands from the optical to 
X-rays. The high mass, intermediate redshift, and compact angu- 
lar scale of these systems make them particularly promising targets 
for the detection of WHIM emission. Previous studies of the large- 
scale environment of clusters at similar redshift (e.g. Kodama et 
al. 2001) have either focused on a single object or on systems of 
relatively low mass. Since more massive clusters are connected to 
more, and denser, filaments than less massive systems (Colberg, 
Krughoff, & Connolly 2005, Colberg 2007), a systematic survey 
of the environment of a statistically complete sample of the most 
massive clusters at these redshifts is advantageous in the context of 
the science goals outlined above. 

The structure of this paper is as follows. In §2 we describe 
the selection of our cluster sample. §3 describes the observations 
and data reduction as well as the methods used in our analysis of 
the data. We present and discuss our results in §4. §5 contains our 
conclusions and identifies areas for future work. Throughout the 
paper we assume a cosmology with Qm — 0.3, = 0.7, and 
h = 0.7. 



2 CLUSTER SAMPLE SELECTION 

The galaxy clusters used in this study are a subset of the sys- 
tems detected in the Massive Cluster Survey (MACS), a survey de- 
signed with the aim of finding the distant counterparts to the most 
X-ray luminous nearby clusters (Ebeling, Edge, & Henry 2001). 
We here focus on the complete sample of the twelve most dis- 
tant MACS clusters at z > 0.5 as published by Ebeling et al. 
(2007). As the most extreme representatives of a sample of al- 
ready extraordinarily massive clusters (typically 1O 15 M ), these 
most distant MACS clusters have been high-priority targets of an 
ongoing comprehensive multiwavelength follow-up campaign, en- 
suring that wide-field optical and X-ray data are available for all 
12 clusters. At the sample's median redshift of z = 0.55 the ex- 
pected scale of the large-scale structure we aim to study is ide- 
ally matched to the field of view of modern wide-field cameras, 
eliminating the need for mosaicking. A pilot study of one of these 
clusters (MACS J0717.5+3745) unveiled a - 4 Mpc long filament 
(Ebeling, Barrett, & Donovan 2004), a discovery that provided the 
impetus for this comprehensive study of the environment surround- 
ing similarly massive clusters in order to determine whether such 
filaments are rare for these systems or whether they are the norm. 



3 DATA ACQUISITION AND ANALYSIS 
3.1 Optical 

As part of our comprehensive observational follow-up campaign of 
MACS clusters we imaged all 12 clusters of this sample with the 
SuprimeCam wide-field imager on the Subaru 8.2m telescope on 
Mauna Kea (Miyazaki et al. 2002). Covering a field of view of ap- 
proximately 34 x 27 arcmin 2 (corresponding to 13 x 10 Mpc 2 at 
the median redshift of our cluster targets), SuprimeCam is ideally 
suited to study the large-scale environment of massive clusters. We 
here use data obtained in the V and R broadband filters, taking ad- 
vantage of the fact that these passbands bracket the D4000 A break 
at the redshifts of our clusters and thus allow an efficient selection 
of early- type cluster members using the cluster red sequence. 

The imaging data used here were obtained in several observ- 
ing runs between December 2000 and November 2005 in variable, 
but usually near-photometric, conditions. The seeing ranged from 
0.7 to 1.2 arcsec as measured from the final stacked images. The 
integration times were, typically, 48 minutes in R and 36 min in V. 
All data were reduced employing standard techniques which have, 
however, been adapted and refined specifically for the analysis of 
SuprimeCam data; details of the procedure are given by Donovan 
(2007). 

In order to allow the measurement of robust V — R colors for 
all objects within the field of view we ensured that the V and R 
band images for each cluster had the same effective spatial resolu- 
tion by blurring the image obtained in better conditions until its see- 
ing, as measured for several hundred stars in the field, matched the 
one obtained in poorer seeing. The seeing in these seeing-matched 
images ranged from 0.9 to 1.2 arcsec for the 12 clusters of our sam- 
ple. For the clusters falling within the footprint of the Sloan Digital 
Sky Survey (SDSS), the photometric zeropoints were derived us- 
ing the stars in the SDSS catalogu^E For those cluster fields not 
covered by the SDSS we obtained cross-calibration in the form of 
3 -sec. exposures on nearby SDSS fields. 

Object catalogues for each cluster were then created with SEx- 
tractor Version 2.4.3 (Bertin & Arnouts 1996) in "dual image" 
mode with the R-band image as the reference detection image. We 
separated stars from galaxies using fits to the stellar sequence in 
both the magnitude vs. peak surface-brightness distribution and the 
magnitude vs. half-light radius distribution. We fit a straight line to 
the stellar locus and consider objects farther than 3a from the best 
fit to be galaxies. For the galaxies thus selected we created V — R 
vs R color-magnitude diagrams using fixed apertures of 1 arcsec 
radius for the color measurement, combined with total isophotal 
magnitudes in the R band as returned by SExtractor's MAG_AUTO 
parameter. 

We then co- added the color-magnitude diagrams for all of the 
clusters after shifting the various data sets along the V — R axis 
such that, for each cluster, the median color of the fifty brightest 
cluster galaxies matched the respective value for a reference clus- 
ter near the median redshift of our sample (MACS J0717. 5+3745, 
z = 0.545). The R magnitudes were corrected to the same me- 
dian cluster redshift of z = 0.55 using the distance modulus and 
a K-correction obtained with a template appropriate for a red se- 
quence galaxy and generated using a suitable model from Bruzual 
& Chariot (2003) with a Salpeter initial mass function, solar metal- 
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system using the equation of Lupton (2005) on the SDSS homepage: 
http://www.sdss.org/dr5/algorithms/sdssUBVRITransform.html 



Large-Scale Structure Around Distant MACS Clusters 3 



2joi 1 1 1 1 1 1 1 1 1 1 1 1 r ~; — r 




IS 20 22 24 

B 



Figure 1. Stacked color-magnitude diagram for our sample, after correcting 
magnitudes to a nominal cluster redshift of z = 0.55 and shifting V — R 
colors such that the median color of the first to fortyninth-ranked galaxies 
in all cluster red sequences match. The orange polygon marks the region 
within which we consider galaxies to be cluster members. 

licity, and assuming a formation redshift of z = 5.0. The re- 
sulting stacked color-magnitude diagram is shown in Fig. 1. The 
color- selected galaxy sample used in our analysis was obtained by 
performing a linear fit to the red sequence in this stacked color- 
magnitude diagram down to an R magnitude of 22.5. We then se- 
lected all galaxies within 2a (as defined as the gaussian width of 
the galaxies around the best-fit line) of this best-fit red sequence 
down to an R magnitude of 23.5 as shown by the orange polygon 
in Fig. 1. This magnitude limit was chosen by measuring the global 
90% completeness limit for all galaxies in the entire samplqj and 
corresponds to 0.03L* at z = 0.55. For the purposes of this work 
we consider all of these red- sequence selected galaxies to be cluster 
members and tracers of the densest and most highly evolved part of 
the galaxy population of our cluster targets and their environments. 
Further details for each cluster are given in Table 1 . 

3.2 X-ray Observations 

The X-ray data were obtained with the ACIS-I instrument (the 
imager of the Advanced CCD Imaging Spectrometer) aboard the 
Chandra X-Ray Observatory between January 2001 and December 
2005. The exposure time ranged from 16 to 70 ks. All observations 
were conducted in VFaint mode and used chips 1, 2, 3, 4, and 6 of 
the ACIS CCD array. The data were reduced using standard pro- 
cedures of the CIAO 3.0.2 software suite. The X-ray images in the 
0.5-7 keV band were smoothed using the adaptive-kernel smooth- 
ing algorithm, ASMOOTH (Ebeling, White, & Rangarajan 2006) to 
produce the contour maps discussed in the next section. 



4 RESULTS AND DISCUSSION 

Figure 2 shows, for each of the twelve clusters in our sample, con- 
tours of the adaptively smoothed, projected surface density of the 

2 For the most distant cluster of our sample, MACS J0744. 8+3927, 90% 
completeness is already reached at R ~ 22 (after correcting to z = 0.55); 
we hence applied this more stringent cut for this one cluster. 



red-sequence selected galaxies as obtained with ASMOOTH. By de- 
sign, all features within the bold outline are significant at the 3a 
confidence level. Figure 3 shows the same data in a scaled rep- 
resentation, but only for the cluster cores and with X-ray surface 
brightness countours overlaid. 

4.1 Optical Morphology 

All but three of the twelve clusters show evidence for at least one 
satellite cluster, the exceptions being MACS J091 1.2+1746 
(which, however, consists of two close subclusters), 
MACS J1423. 8+2404, and MACS J2129.4— 0741. The promi- 
nent satellite cluster almost 10 arcmin southwest of the core 
of MACSJ0018+1626 (A.K.A. C10016+16, Spinrad 1980) has 
been previously identified (Connolly et al. 1996) but the filament 
connecting the two systems is reported here for the first time. More 
than half of all systems feature more than one apparent satellite 
(MACS J0018+1626, MACS J0025.4-1222, MACS J0257.1- 
2335, MACS J0454. 1-0300, MACS J0647.7+7015, 

MACS J0717.5+3745, MACS J2214.9- 1359) and 

an equal number show well-defined filamen- 
tary structure (MACS J00 18+ 1626, MACSJ0257.1- 
2335, MACS J0454. 1-0300, MACS J0647.7+7015, 

MACS J0717.5+3745, MACS J2129.4-0741, 

MACS J2214.9- 1359). The filament extending to the south- 
east of the main cluster in MACS J07 17.5+4735 was previously 
identified by Ebeling et al. (2004) who point out that its orientation 
matches that of the merger axis of the cluster proper, suggesting 
that the system has grown for many Gyr via infall of matter 
along this filament. We find similar alignments between filaments 
and the galaxy-density contours of the cluster core for almost 
all other systems exhibiting filamentary structure (seven out of 
twelve as listed above), providing strong observational evidence 
that the growth of most clusters occurs along a single preferred 
direction, be it through major mergers and/or steady accretion, and 
that these "assembly lines" are stable over remarkably long time 
scales. Even one cluster without prominent large-scale filaments 
(MACS J0025.4- 1222) exhibits an elongation of the outermost 
contours of the galaxy density near the main cluster that appears to 
point toward a satellite cluster, suggesting that an active filament 
may have been present in its direction at earlier times. 

The most complex systems in terms of large-scale structure in 
this sample are MACS J0257. 1-2335, MACS J0647.7+7015, and 
MACS J2214.9-1359. Each of these three features at least three 
collapsed cluster cores separated by several Mpc, as well as two 
nearly orthogonal axes of merging and infall along filaments. The 
large-scale morphology of the galaxy distribution in these fields 
thus strongly suggests that the respective MACS clusters occupy 
intersections of multiple filaments in the cosmic web. The detection 
of several filaments connecting satellite clusters to the main cluster 
provides strong observational support for the results of numerical 
simulations (Colberg et al. 2005; Colberg 2007) which suggest (a) 
that cluster pairs with very close separations (< 5h _1 Mpc) are al- 
ways connected by filaments, and (b) that, statistically, the number 
of filaments increases with cluster mass. 

A closer look at the prevalence of substructure in the cluster 
cores - indicative of an ongoing merger - and the strength of the 
filamentary extensions connecting these cores to their large-scale 
surroundings provides a fascinating snapshot of the wide variety 
in the dynamics and geometry of cluster growth. MACS J2214.9- 
1359, for instance, exhibits only subtle signs of disturbance in its 
core region (primarily a 90-degree rotation in the mild ellipticity 
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Table 1. MACS Clusters at z > 0.5 



Cluster Name a 


Redshift 


# of Galaxies 6 


*exp (V, R) (s) c 


Date Observed (V, R) 


t ex p (X-ray) (ks) c 


MACS J0018.5+1626 


0.5476 


1515 


2160, 2880 


Sep. 03, Dec. 02 


70.1 


MACS J0025.4-1222 


0.5840 


1170 


2160, 2880 


Sep. 03, Dec. 02 


27.4 


MACS J0257.1-2335 


0.5058 


883 


2160, 5280 


Dec. 02, Dec. 00 


20.5 


MACS J0454. 1-0300 


0.5377 


997 


2160, 3240 


Nov. 05, Mar. 05 


16.9 


MACS J0647.7+7015 


0.5912 


1478 


2160, 2880 


Sep. 03, Feb. 04 


22.1 


MACS J0717.5+3745 


0.5457 


1872 


2160, 2880 


Dec. 02, Dec. 00 


61.7 


MACS J0744.8+3927 


0.6984 


1231 


1440, 2880 


Apr. 03, Dec. 02 


53.2 


MACS J0911.2+1746 


0.5041 


849 


2160, 2880 


Apr. 03, Apr. 03 


26.7 


MACS Jl 149.5+2223 


0.5445 


1589 


2160, 2880 


Apr. 03, Apr. 03 


22.0 


MACSJ1423.8+2404 


0.5432 


1093 


2160, 2400 


Jun. 02, Jun. 02 


20.2 


MACS J2129.4-0741 


0.5347 


1108 


1440, 2880 


Jun. 02, Jun. 01 


21.7 


MACS J2214.9-1359 


0.5030 


1462 


2160, 2880 


Jul. 07, Jul. 03 


20.1 



a For coordinates and other cluster properties, see Ebeling et al. 2007 

b Number of galaxies selected along the red sequence using the method described in section 3.1 
c Exposure times of the optical and X-ray data in seconds and kiloseconds, respectively 



of the galaxy-density isophotes from the very core to the outer re- 
gions) suggesting a state of advanced relaxation since the last ma- 
jor merger. On larger scales, however, faint filaments still trace the 
paths of past and future infall of matter. By stark contrast, the dis- 
turbed core of MACS J0647.7+7015 is clearly the scene of ongo- 
ing merging, and the strength of the filemantary connection to a 
series of satellite clusters to the southeast of the main cluster sug- 
gest that this phase of violent growth is to continue for many Gyr to 
come. Viewed as a whole, our twelve snapshots of the large-scale 
structure around the most massive galaxy clusters at z ~ 0.55 thus 
corroborate work by Colberg et al. (1999) whose numerical simula- 
tions show strong correlations between infall patterns and the large- 
scale structure surrounding massive clusters, with preferred direc- 
tions for the assembly of clusters being established by tidal forces. 
Our observational results strongly suggest that, although multiple 
filaments are often present, one of them usually dominates in fun- 
neling matter onto the cluster occupying the respective vertex of 
the cosmic web. 

Finally, a minority of our cluster targets exhibit no- 
ticeably elongated galaxy-density contours in their core 
regions (MACS J0744.8+3927, MACS Jl 149.5+2223, 
MACS J1423. 8+2404), but little sign of substantial filaments 
on larger scales. Since the orientation of major filaments is random 
with respect to our line of sight it is possible that significant 
large-scale structure is present but not visible in projection. Indeed 
extensive spectroscopic follow-up of the galaxy population of 
these systems finds strong evidence for line-of- sight substructure 
in one of these three, MACS Jl 149.5+2223 (Smith et al. 2008, in 
preparation). At least one cluster in our sample can, however, be 
considered to be near virialization and safe from major mergers for 
some time, namely MACS J 1423. 8+2404 (see also the following 
section). 



4.2 X-ray Emission 

Although the small field of view, compared to SuprimeCam, of 
ACIS-I allows us to study only the X-ray emission from the main 
clusters shown in Fig. 2, fresh insights about the dynamical state 
of these clusters can be obtained by comparing the X-ray surface- 
brightness contours to the optical morphology of the main cluster 



(Fig. 3). In all cases we find the orientation of any elongation in the 
X-ray contours around the cluster core to be well aligned with the 
direction of infall as traced by the large-scale galaxy distribution. 
This is expected since the intra-cluster gas responds to the same 
gravitational forces as the galaxies and should thus, globally, show 
similar bulk motions. 

There is a fundamental difference between gas and galax- 
ies, however, in that the former is a strongly interacting viscous 
medium, whereas the latter are, on cluster scales, weakly or non- 
interacting point- like objects. As a consequence, gas and galaxies 
can show noticeable segregations in non-equilibrium situations like 
major cluster mergers. Indeed we find two systems for which the 
centroids of the gas and the galaxy distributions in the cluster cores 
are noticeably displaced from each other, namely MACS J0025.4- 
1222 and MACS J0717.5+3745, with the displacement correspond- 
ing to roughly 200-300 kpc in the plane of the sky. In either case the 
X-ray emission peaks between two well separated galaxy concen- 
trations. A third, less obvious, example of the same phenomenon 
is MACS J091 1.2+1746, a cluster consisting of two components of 
very different mass about 1 Mpc apart. A comparison of the galaxy 
and gas distributions around the less massive component shows the 
centroid of the X-ray emission to be significantly displaced from 
the peak in the galaxy surface density. The cause is, most likely, the 
same for all three systems mentioned, namely an ongoing merger in 
which the galaxy populations of the two merging clusters have al- 
ready passed through each other with little interaction, whereas the 
collision of the dense gas in the cluster cores creates shocks through 
which kinetic energy from the bulk motion of the gas is partly con- 
verted into thermal energy and partly redirected into transverse ve- 
locity components. 

Similar cases of a clear segregation between gas and galax- 
ies have been observed in other clusters and have helped greatly to 
improve our understanding of cluster mergers. A particularly spec- 
tacular example is the "bullet cluster" 1E0657-558 (Markevitch 
et al. 2002) in which the fortuitous alignment of the merger axis 
with the plane of the sky enabled remarkably detailed and quantita- 
tive studies of the dynamics of the merging process. It remains, so 
far, the only system to allow direct constraints to be placed on the 
dark matter self-interaction cross section (Markevitch et al. 2004, 
Clowe et al. 2006). The potential implications of significant off- 
sets between the centroids of the gas, galaxy, and dark-matter dis- 
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tribution in clusters are thus far-reaching, and further follow-up 
studies, in particular a gravitational-lensing analysis of the total- 
mass distribution, of MACS J0025.4-1222, MACS J07 17.5+3745, 
and MACS J091 1.2+1746 may offer a rare opportunity to test the 
results obtained for 1E0657-558 on the nature and fundamental 
properties of dark matter. 

Two of our 12 MACS clusters (MACS J0744. 8+3927 and 
MACS J1423. 8+2404) deserve special mention as the only 
ones showing obvious cooling cores in our X-ray images. 
MACS J1423. 8+2404 is also the sole cluster that shows neither 
obvious large-scale structure nor significant substructure in the 
core region. Although, again, structure along the line of sight can- 
not be ruled out conclusively, the system's optical and X-ray ap- 
pearance suggests that it is close to virialization (see also Ebel- 
ing et al. 2007). While a cooling core is also clearly present in 
MACS J0744. 8+3927, the large-scale environment of this system 
and thus its overall assembly history and future are, at present, dif- 
ficult to assess owing to the limited depth of the existing Suprime- 
Cam data for this most distant cluster in our sample. 



5 SUMMARY 

We use the projected surface density of galaxies near the cluster 
red sequence to trace the large-scale environment around all twelve 
MACS clusters at z > 0.5 from the sample of Ebeling et al. (2007). 
Since the galaxy-color criterion applied selects primarily early-type 
galaxies dominated by old stellar populations this approach is well 
suited to map the galaxy distribution in regions of moderate to high 
density (cluster cores, satellite clusters, and the densest regions of 
the filamentary structures connecting them), but will be less effi- 
cient in unveiling low-density structures such as galaxy sheets or 
the environment customarily referred to as the "infall region", i.e., 
the region around the cluster-field interface. Our study uses wide- 
field imaging data obtained with the SuprimeCam imager on the 
Subaru 8m telesope on Mauna Kea and probes scales of typically 
10 Mpc at the median redshift of z = 0.55 of our cluster targets. 
By correcting the distribution of observed galaxy magnitudes to a 
common reference redshift (again z = 0.55) and applying a global 
magnitude cut we ensure that all systems are sampled uniformly to 
0.03L^, with the exception of MACS J0744.8+3927 for which our 
existing data become incomplete at 0.2L^. 

We find the vast majority of these most massive clusters at 
z ~ 0.55 to be embedded in a highly complex large-scale environ- 
ment characterized by satellite groups/clusters and extensive fila- 
ments. In all cases the alignment of any substructure and/or ellip- 
ticity of the galaxy distribution in the cluster cores with the fila- 
ments found on much larger scales strongly suggests that the ongo- 
ing growth of these systems occurs through merging and infall of 
matter along a single preferred axis that is stable over many Gyr. 
Several clusters are found to reside at the intersection of multiple 
distinct axes of infall which appear to have actively contributed to 
the assembly of the cluster in the past, presumably alternating in 
dominance. The observed large-scale morphologies and their im- 
plications for the dynamics of cluster formation and evolution thus 
provide strong observational support for the results of numerical 
work that finds the growth of the most massive galaxy clusters to 
occur at the densest vertices of the cosmic web, tightly linked to 
filaments funneling matter onto the central cluster. Since this work, 
by design, is biased toward detecting the denser regions of the en- 
vironments studied, the extensive filaments found should be well 
suited for observational attempts to directly detect and study the 



properties of the "Warm-Hot Intergalactic Medium" (WHIM) hy- 
pothesized to contain a significant fraction of the baryons in the 
universe. 

Chandra observations of all clusters in our sample allow a 
comparison of the distribution of galaxies with that of the gas in the 
central cluster regions. We find the morphology of gas and galax- 
ies to be similar on large scales, consistent with the expectation 
that in near-equilibrium conditions all cluster components partic- 
ipate similarly in gravitation-induced bulk motions. We, however, 
find significant offsets between the centroids of the two distribu- 
tions in the cases of MACS J0025.4-1222, MACS J07 17.5+3745, 
and MACS J091 1.2+1746. All of these systems are thus promis- 
ing targets for in-depth follow-up studies of the distributions of 
gas, galaxies, and dark matter in merging clusters that hold great 
promise for advancing our understanding not just of the physics 
and dynamics of cluster mergers but specifically also of the nature 
and properties of dark matter (Bradac et al. in preparation). 

Finally, the wide range of environments discovered here 
confirms MACS clusters as ideal targets for detailed studies of 
the impact of environment on galaxy evolution. First results on 
the distribution and origin of E+A galaxies from an extensive 
spectroscopic follow-up campaign of galaxies in the field of 
MACS J07 17.5+3745 underline the significant advantage of being 
able to sample the full range of environments, from the low-density 
field to the extremely dense cluster cores, in a single observational 
target (Ma et al. 2008). 

We would like to thank Leif Wilden and Glenn Morris of Stan- 
ford University for assistance with the reduction of some of the 
optical data used in this study. HE gratefully acknowledges finan- 
cial support from NASA LTSA grant NAG 5-8253 and SAO grants 
G02-3168X, G03-4168X, and G03-4164X. 



REFERENCES 

Bertin E., Arnouts S. 1996, A & AS, 117, 393 
Bruzual G., Chariot S. 2003, MNRAS, 344, 1000 
Cen R., Ostriker J. P. 1999, 514, 1 
Clowe D. et al. 2006, ApJL, 648, L109 
Colberg J. M. et al. 2000, MNRAS, 319, 209 
Colberg J. M., Krughoff K. S., Connolly A. J. 2005, MNRAS, 
359, 272 

Colberg J.M., White S.D.M., Jenkins A., Pearce F.R. 1999, MN- 
RAS, 308, 593 
Colberg J.M. 2007, MNRAS, 375, 337 

Connolly A.J., Szalay A. S., Koo D., Romer A. K., Holden B., 

Nichol R.C., Miyaji R. 1996, ApJL, 473, L67 
Dave R., Hernquist L., Katz N., Weinberg D. H. 1999, ApJ, 511, 

521 

Dave R. et al. 2001, ApJ, 552, 473 
Donovan D. 2007, PhD Thesis, University of Hawaii 
Ebeling H., Barrett E., Donovan D. 2004, ApJL, 609, L49 
Ebeling H., Edge A. C, Henry J. P. 2001, ApJ, 471, 694 
Ebeling H., White D.A., Rangarajan F.V.N 2006, MNRAS, 368, 
65 

Ebeling H. et al. 2007, ApJL, 661, 33 
Einasto J. et al. 2007, A&A, 462, 811 
Erdogdu P. et al. 2004, MNRAS, 352, 939 
Evrard A. E., et al. 2002, ApJ, 573, 7 
Hallman E. J. et al. 2007, ApJ, 671, 27 
Kodama T., Arimoto N. 1997, A&A, 320, 41 



6 S. Kartaltepe et al. 



Ma C.-J., Ebeling H., Donovan D., Barrett E. 2008, ApJ, submit- 
ted 

Markevitch M., Gonzalez A.H., David L., Vikhlinin A., Murray 
S., Forman W., Jones C, Tucker W. 2002, ApJ, 576, L27 

Markevitch M., Gonzalez A.H., Clowe D., Vikhlinin A., Forman 
W., Jones C, Murray S., Tucker W. 2004, ApJ, 606, 819 

Miyazaki S. et al. 2002, PASJ, 54, 833 

Nicastro F, Elvis M., Fiore F, Mathur S. 2005, AdSpR, 36, 721 
Peacock J.A. et al. 2001, Nature, 410, 169 
Soltan A.M., Freyberg M.J., Hasinger G. 2005. A&A, 436, 67 
Spinrad H. 1980, IAUS, 92, 39 

Zandivarez A., Merchan M.E., Padilla N.D. 2003, MNRAS, 344, 
247 



Large-Scale Structure Around Distant MACS Clusters 7 




19 m 00 s 30 s 18 m 00 8 

Right Ascension (2000) 




26 m 00 s 30 s 25 m 00 s 

Right Ascension (2000) 



MACSJ0257 




58 m 00 s 


30 s 57 m 00 s 30 s 
Right Ascension (2000) 


02 h 56 m 00 s 







r 


MACSJ0647 






















- 


1 1 





MACSJ0454 





54 m 00 s 
Right Ascension (2000) 



04 h 53 m 00 s 




50 49 



Right Ascension (2000) 



47 46 06 h 45 m 



30 s 18 m 00 s 30 s 17 m 00 s 07 h 16 m 30 s 

Right Ascension (2000) 



Figure 2. Contours of the the adaptively smoothed projected surface-density distribution of galaxies on the cluster red sequence. The contours are logarithmi- 
cally spaced by 20%. Features within the bold line are significant at the 3a confidence level. 
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Figure 2 - continued 
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Figure 3. Scaled representation of the same data as in Fig. 2 but showing only the core regions of each cluster. Overlaid in blue are the contours, logarithmically 
spaced by 20%, of the adaptively smoothed X-ray emission from the diffuse intra-cluster medium as observed with Chandra/ACIS-I in the 0.5-7 keV band. 
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Figure 3 - continued 



